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Abstract 

This paper deals with stability of analytical and numerical solutions of nonlinear stochastic 
delay differential equations (SDDEs). Sufficient conditions for stability in mean square, 
contractivity in mean square and asymptotic contractivity in mean square of the solutions 
of nonlinear SDDEs are given, which provide a unified theoretical treatment for constant 
delay and variable delay, bounded delay and unbounded delay. Stability in mean square, 
contractivity in mean square and asymptotic contractivity in mean square of the backward 
Euler method are investigated. It is shown that the backward Euler method can inherit the 
properties of the underlying system. The main theorems in this paper are different from 
Razumikhin-typc theorems. It is not required to construct or find an appropriate Lyapunov 
functional when the theorems in this paper are applied. 

AMS subject classification: 34K50, 65C20, 65C30. 
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1 Introduction 

Many real-life phenomena in physics, engineering, economics, etc can be modeled by stochastic 
differential equations (SDEs). The rate of change of such a system depends only on its present 
state and some noisy input. However, in many situations the rate of change of the state depends 
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not only on the present but also on the past states of the system. Stochastic functional differential 
equations (SFDEs) give a mathematical formulation for such system. For more details on SFDEs, 
the reader may refer to \12\ [T3] and the references therein. 

SFDEs also can be regarded as a generalization of deterministic functional differential equa- 
tions if stochastic effects are taken into account. For deterministic Volterra functional differential 
equations (VFDEs), Li [3] discussed stability, contractivity and asymptotic stability of the solu- 
tions in Banach spaces. The main techniques used in [8] are as follows: The author introduced 
a so-called ^-perturbed problem and constructed an auxiliary function Q(t). The former is 
beneficial to deal with a wide variety of delay arguments and the latter is the key to obtain 
main results in that paper. [8] provided a unified framework for stability analysis of solutions to 
nonlinear stiff problems in ordinary differential equations (ODEs), delay differential equations 
(DDEs), integro-differential equations (IDEs) and VFDEs of other type appeared in practice. 
In 2008, Wang [T7] extended successfully the theory in [8] to nonlinear Volterra neutral func- 
tional differential equations (VNFDEs). Wang and Zhang [18] proved implicit Euler method 
can preserve the stability of VFDEs and VNFDEs. 

It is natural to ask whether the solutions of SFDEs possess the properties similar to those 
presented in [S] and which methods can reproduce the properties. Due to features of the stochas- 
tic calculus the numerical analysis of SFDEs differs in some key areas from the already well- 
developed area of the numerical analysis of their deterministic counterparts. In the literature, 
much attention on numerical stability has been focused on a special class of SFDEs, that is, 
stochastic delay differential equations (SDDEs). The stability includes mean-square stability, 
asymptotic stability, exponential mean-square stability and so on. Some results on numerical 
stability can be found in [T| 110 1 H9 | f20] , 121 j . However, most of the existing numerical stability 
theory of SDDEs only deal with bounded lags. Far less is known for long-run behavior of SDDEs 
with unbounded lags. Very recently, Fan, Song and Liu [3] discussed the mean square stability 
of semi-implicit Euler methods for linear stochastic pantograph equation. Moreover, to the best 
knowledge of authors, there is no work on contractivity analysis of numerical methods for SD- 
DEs. These motivate our work. This is the first paper to investigate stability and contractivity 
of general SDDEs with bounded and unbounded lags and study whether numerical methods can 
preserve the properties. 

The main contributions of this paper could be summarized as follows: 

(a) Sufficient conditions for stability in mean square, contractivity in mean square and asymp- 
totic contractivity in mean square of the solutions of nonlinear SDDEs are given, which provide 
a unified theoretical treatment for constant and variable delay, bounded and unbounded delay. 
Applicability of the theory is illustrated by linear and nonlinear SDDEs with a wide variety 
of delay arguments such as constant delays, piecewise constant arguments, proportional delays 
and so on. The theorems established in this paper work for some SDDEs to which the existing 
results cannot be applied. The main results of analytic solution in this paper can be regarded as 
a generalization of those in [5] restricted in finite-dimensional Hilbert spaces and finitely many 
delays to the stochastic version. 

(b) It is proved that the backward Euler method can reproduce stability in mean square, 
contractivity in mean square and asymptotic contractivity in mean square of the underlying 
system. In particularly, Theorem 14.21 and Theorem 14.41 show that the backward Euler method 
can inherit the contractivity and the asymptotic contractivity without any constraint on stepsize. 

We point out that the main theorems in the present paper are different from the Razumikhin- 
type theorems in the literature [TJ [11] . It is not required to construct or find an appropriate 
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Lyapunov functional when the theorems in this paper are applied to prove the stability, in 
contrast to the Razumikhin-type theorems. In this sense, it is more convenient to use the 
theorems in this paper than the Razumikhin-type theorems. 

The rest of the paper is arranged as follows. In section 2, we introduce some necessary nota- 
tions and assumptions. In section 3, some criteria for the stability in mean square, contractivity 
in mean square and asymptotic contractivity in mean square of solutions of nonlinear SDDEs 
are established. The main results obtained in this section are applied to SDDEs with bounded 
and unbounded lags, respectively. In section 4, sufficient conditions on stability in mean square, 
contractivity in mean square and asymptotic contractivity in mean square of the backward Euler 
method are derived. Stability of analytic and numerical solutions of SDDEs with several delays 
is discussed in section 5. 



2 Stochastic delay differential equations 

Let (0, J 7 , {J^j^a, P) be a complete probability space with a filtration {Ft}t>a satisfying the 
usual conditions (i.e., it is right continuous and T a contains all the P-null sets). Let w(t) := 
(wi(t), ...,w m (t)) T be an m-dimensional Wiener process defined on the probability space. And 
| • | is used to denote both the norm in complex Hilbert space C d and the trace norm (F-norm) in 
C dxm . Also, C([ti, t 2 \, C d ) is used to represent the family of continuous mappings i/j from [t\, t 2 ] 
to C d . Finally, L p -p ([ii, £2]; C d ) is used to denote a family of Tt -measurable, C([t\, t 2 \; C d )-valued 
random variables ip = {i[>(u) : t\ < u < £2} such that [IV'He := su Pti<«<t 2 

E\i/)(u)\p < 00. E 

denote mathematical expectation with respect to P. 

Consider the following initial value problems of SDDEs in the sense of Ito 

dx(t) = f(t, x(t),x(t - T{t)))dt + g(t, x(t),x(t - r(t)))dw(t), t £ [a, b], (2.1a) 
x(t)=£(t), te[a-r,a], (2.1b) 

where a,b,T are constants, —00 < a < b < +00, r > 0, r(t) > 0, inf (t — r(t)) > a — r, 

a<t<b 

f S L P Ta {[a - T,a];C d ),p > 2, / : [a, b] x C d x C d -»■ C d ,g : [a,b] x C d x C d -> C dxrn are given 
continuous mappings. We assume that / and g satisfy the following conditions. 

For each R > there exists a constant Cr, depending only on R, such that 

(2.2) 

\f(t,xi,y) - f(t,x 2 ,y)\ < C R \x\ - x 2 \, V \x 2 \ V \y\ < R, 

SR(a;i - x 2 ,f(t,xi,y) - f(t,x 2 ,y)) < a(t)\xi - x 2 \ 2 , (2.3) 
\f(t,x, yi ) - f(t,x,y 2 )\ < 0(t)\yi - y 2 \, (2.4) 
\g(t,xi,yx) - g(t,x 2 ,y 2 )\ < -yi(i)[cci - x 2 \ +72(^)^1 — 2/2 1 , (2.5) 

for all t E [a,b],x,xi,x 2 ,y,yi,y 2 G C d . Here a(t), (3(t), Ji(t) and 72 (t) are continuous real- valued 
functions. We introduce the following notations: 

/4 0) = inf r(t)>0, 4 0) (ti,t 2 )= mf (t - r(t)) > a - r, 

a<t<b ti<t<t2 

Vti,t 2 : a < t\ < t 2 < b. 

For convenience, we denote by ST>(a, (3, 71 , 72) the all problems (|2.ip which satisfy the conditions 
(|2.2|) - (|2.5p . Examples of linear and nonlinear equations which satisfy the conditions (|2.2p - (|2.5|) 
are given in the next section (see Example 13.161 and Example I3.20p . 
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In order to deal with a wide variety of delay arguments in this paper, we introduced the 
so-called ^-perturbed problem of (|2.1|) which first introduced by Li[8] for VFDEs. Let 

f n (t,x,i/>) = /(t,x,V {n,i) ), 9n(t,x,ilt) = g(t,x,^ n '% (2.6) 
for t E [a,b],x £ C d , ip € C([a - r, b];C d ), where V (n,i) is defined by 

*<*<•>-{ «?-i>. :il::i',V"'' (2 - 7) 

Then the initial problem 

dx(t) = f n (t,x(t),x(t - r(t)))dt + g n (t,x(t),x{t - T(t)))dw(t), t e [a,b], (2.8a) 
x(t) =£(t), te[a-T,a] (2.8b) 

is said to be the —perturbed problem of the problem (|2.ip . Here the natural number n > ~ 
can be arbitrarily given. Notice that for the case of r = 0, we can replace it by some positive 
number f and define £(u) = £(a) for u € [a — f,a\. Therefore, without lose of generality, we 
always assume r > . 

It is easy to verify that, if problem (|2.ip £ ST>{a, f3, 71, 72), then its —perturbed problem 
USSJ € SV(a,P,~/i,~/ 2 )- Let f(t) = max{r(t), £}, we have 

/iS 0) = inf f(t) > i. (2.9) 

a<t<b n 

3 Stability analysis of SDDEs 

We will discuss three types of stability of SDDEs in this section and here give definitions of 
them. 

Definition 3.1 The solution of problem 112. 1\) is said to be stable in mean square if 

nx(t) - y(t)\ 2 < C sup E\m-r}(0)\ 2 , (3.1) 

a~r<9<a 

where y(t) is a solution of any given perturbed problem 

f dy(t) = f(t, y(t),y(t - r(t)))dt + g(t, y(t),y(t - r(t)))dw(t), t G [a, b], (3.2a) 
\y(t) = V (t),te[a-r,a], 77 € i£ s ([a - r, a]; C d ). (3.2b) 



Definition 3.2 T/ie solution of problem \2. 1\) is said to be contractive in mean square if US. 1\) 
with C < 1 holds. 

Definition 3.3 The solution of problem \2.1\) is said to be asymptotically contractive in mean 
square if 



lim E\x(t) - y(t)\ z = 0, 

t— S-+QO 



where [a, b] replaced by [a, +00) 
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Remark 3.4 In the strict sense, the solution of 12.1]) should be called generalized contractive in 
mean square if \3.1\j with C < 1 holds. For convenience, we call it contractive in mean square 
for brevity. 

There exist well-known stability definitions in the literature which are closely related to those 
presented in this paper, but there are differences among them. The existing notions of stability 
include mean-square stability for stochastic differential equations, that is, lim E|x(t)| 2 =0 (cf. 

t— >+oo 

114V i exponential mean-square contraction of trajectories for stochastic differential equations 
with jumps (cf. f^). The contractivity in mean square is weaker than that in f^. 

The continuity of / and g implies that 

|/(i,0,0)| <C, \g(t,0,0)\ <C, t€[a,b], (3.3) 

where C only depends on f,g and the interval [a, b]. We note that condition (|2.5p implies that 
the diffusion coefficient g satisfies local Lipschitz condition 

\g(t,xi,yx) - g(t, x 2 ,y 2 )| < C R (\xi - x 2 \ + \yi - 2/2I) , 

t G [a, b],xi,x 2 ,yi,y2 G C d , \xi\ V \x 2 \ V \yt\ V \y 2 \ < R. 

In fact, we can choose any Cr with Cr > max{ max 71 (t), max 72 (t)}- Using (|2.3|) . (|2.4|) and 

a<t<b a<t<b 

, we have 

f(t, x, y)) = K(x - 0, fit, x, y) - f(t, 0, y) + f(t, 0, y) - f(t, 0, 0) + f(t, 0, 0)' 

< a(t)\x\ 2 + (3{t)\x\\y\ + C\x\ < Ci + |x| 2 + |y| 2 ) , (3.5) 

where C\ only depends on C, max a(t) and max /3(i). By ()2.5|) . we have 

a<t<b a<t<b 

\g(t,x,y)\ 2 < 2\g(t,x,y) - g(t, 0,0)\ 2 + 2\g(t,0,0)\ 2 <C 1 (l + \x\ 2 + |y| 2 ), (3.6) 
where C% only depends on C, max 7 2 (t) and max 7|(t)- 

a<t<b a<t<b 

3.1 Finite interval 

In order to prove the main theorems in this section, we prepare the following lemmas. 

Lemma 3.5 Assume that problem 112. 1\) G SV(a, [3, 71, 72). Then for each p > 2 there is 
= C(p, a, b, a, [3, 71,72) such that 

E( sup |x(t)| p ) < (7(1 +E( sup =: A. (3.7) 

a<t<b a—T<t<a 

Proof. For every integer k > 1, define the stopping time 

p k = inf{* G [a, 6] : sup |x(6>)| > fc}, (3.8) 

a-r<6»<t 
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where we use the convention inf = b. Clearly, py. f b almost surely as k — >■ +00. Let = 
x(t A pi~). Using the ltd formula, we have for a <t <b 



l + [^)| 2 ) § = (l + |£(a)^ f 



+ P j (l + |^( S )| 2 ) 2 M(x k (s)J(s,X k (s),X k ( S -T(s))))I M ds 

+| f (l+|x fc ( S )| 2 ) V | ff ( S ,x fe ( S ),x fe ( S -r( S )))| 2 / M d S 

+ p(p_2) ^* (l + l^^p^^j^fc^))^^ ^( s ),^( s -^( s )))|2 T[[ap&]]rfs 

ft p— 2 

+p/ (l + \x k (s)\ 2 ) 2 R(x k (s),g(s,x k (s),x k (s-T(s))))I [[a , Pk]] dw(s). 



It follows from (|3.5p . ()3.6|) and this expression that 



l + |x fc (t)| 2 ) 2 < 2^(1 + 1^)1^+2^1/" (l + \x k (s)\ 2 ) 2 fl+ sup |x fe (u)| 2 )(fe 

y V y V y V s-r(s)<«<s J 

+pd f (l + |x fc (s)| 2 ) V (l+ sup |x fc (u)| 2 )ds 

ia ^ ' ^ s-t(s)<u<s ' 

ft p-4 

+p(p-2)C 1 (l + \x k (s)\ 2 ) 2 \x k (s)\ 2 (l+ sup |x fe (u)| 2 )ds 



Erg 



s— t(s)<u<s 

+pj (l + |x fc ( S )| 2 ) 2 K(x k (s),g(s,x k (s),x k (s-T(s))))I M dw(s) 
< 2^fl + |e(a)| p ) +Kp+l)Ci / (l+ sup |x fc (u)| 2 ) 2 <fe 



which yields 



E sup (1 + [x fc (s)| 2 )2 < 2^ (l +E sup |£(s)| p ) + C 2 E/ a * (l+ sup |x fc (u)| 2 Vds 

a<s<t ^ a— r<s<a ' ^ a— r<u<s ' 



+pE( sup £ (1 + |x^n)| 2 )" 2 K(x fc (n), 5 ( U ,x fc ( U ),x fc (^-r(n))))/ [Kpfc]] ^(u 

a<s<t 



Applying the Burkholder-Davis-Gundy inequality to the third term at the right hand side of 
this expression, we obtain the bound 



s 



pE[ sup^ / (1 + |x fc (u)| 2 ) P22 »(x fc (n), 5 (n,x A Xn),x fc (n-r(n))))/ [Kpfc]] du;(u) 



< CpE^j (l + \x k (u)\ 2 ) p 2 \x k {u)\ 2 \g{u,x k {u),x k {u-T{u)))\ 2 drf 

< C P E( sup (1 + \x k {u)\ 2 )* f (1 + \x k {u)\ 2 )^\x k (u)\ 2 \g{u,x k (u),x k (u-T{u)))\ 2 d u y 

^a<u<t J a ' 

1 v C 2 — 

< -E sup (l + \x k (u)\ 2 y + ^-2C{E / (l+ sup \x k (s)\ 2 ) 2 du. 
Consequently, 

pt P 

E sup (l + \x k (s)\ 2 ) ^ < 2§(l + E sup |£(s)| p ) + C 3 E / (l + sup \x k {s)\ 2 Ydu 

a<s<t a—r<s<a Ja ^ a—r<s<u ' 

ft p 

= 22(l + E sup |C(s)r) +C 3 E / sup (l + \x k {s)\ 2 Y du. 

a—r<s<a Ja a—T<s<u ^ ' 

Further, we notice that 

E sup (1 + |x fc (s)| 2 ) 2 < E sup (l + |C(s)| 2 )2 +E sup (l + |x fc (s)| 2 ) 2 



a— r<s<t a—r<s<a a<s<t 

< 2 £ 2 1 (1 + E sup |£(s)| p )+E sup (1 + \x k {s)\ 2 y. 

a—r<s<a a<s<t 

Therefore, 

E sup (l + \x k (s)\ 2 y < |2f(l+E sup |^(s)H + C 3 f E sup (l + \x k (s)\ 2 ) * du. 

a—r<s<t £ a—r<s<a Ja a—r<s<u ^ ' 

Now the Gronwall's inequality yields that 

E( sup \x k (s)\ p ) <E( sup (l + |x fc (s)| 2 )i) < |2i(l + E sup |£(s)| p )e C3( *- a) . 

Letting k — > +oo and applying the Fatou's lemma, we obtain the desired result. 

Lemma 3.6 Assume that problem 112.1)) € ST>(ot, f3, ji, 72). TTien i/tere exists a unique solution 
x(t) to equation 12.1)) . 



Proof. Using (j2.2|) . (|2.4[) . (|3.4p . ()3.5p and (|3.6p . we can complete the proof in a similar manner 
to those of Theorem 5.2.7 and Theorem 2.3.5 in 1111. 



Lemma 3.7 Assume that problem 112.1)) € ST>(a, (8,71,72)- TTien we /i<roe 

lim sup E\x(t) -x(s)\ 2 = 0. (3.9) 



t-s-X) 



a<s<t<b 



Proof. The proof uses the techniques employed in the proof of Theorem 2.2 in [5]. Integrating 
(|2.ip gives for a < s < t < b 



(3.10) 



x(t) — x(s) = / f(u,x(u),x(u — T(u)))du+ / g(u,x(u),x(u — r(u)))dw(u). 

J s J s 

Let e(s, t) = x(t) — x(s), 

PH = mf{t G [a, b] : sup \x(9)\ > R}, 

a-T<d<t 

where inf = b. Using the Young inequality: for r _1 + q~ l = 1 

ab < -a r + -^b q , Va,b,6> 
r qdil r 

and letting r = | , q = , we thus have for any 5 > 

E(|e( S ,t)| 2 ) = E{\e(s,t)\ 2 I {pR > b} ) + E{\e(s,t)\ 2 I {pR<b} ) 
< E{\e(s,t)\ 2 I {pR > b} ) + jE{\e{s,t)\*>) + ^^Hpr < b), 



where p > 2. It follows from ()3.7p that 



E(|e(s,t)| p ) < 2 P ~ 1 E( sup \x(s)\ p + sup \x{t)\ p )) < 2 P A. 

a<s<b a<t<b 



We then have 



E(\e(s,t)\ 2 ) <E(\e(s,t)\ 2 I {pR > b} ) + 



2 P+1 5A (p - 2)A 



p p p/(p-2) R p- 

Further, using the Holder's inequality and the ltd isometry we obtain 
E{\e(s,t)\ 2 I {pR > b} ) 



= E 
< 2E 



f(u,x(u),x(u — T(u)))du+ / g(u,x(u),x(u — T(u)))dw(u) 



l {PR>b} 



(3.11) 



(3.12) 
(3.13) 

(3.14) 




f(u,x(u),x(u - T(u)))du\ +\J g(u,x(u),x(u - T(u)))dw(u)\ )I {pR > b} 

< 2(t - s)E^I {pR > b} J \f(u,x(u),x(u-T(u)))\ 2 du) +2EJ \g{u,x(u),x(u - r(u)))\ 2 du 

< 2{t - s)E(l {pR > b} J \f(u,x(u),x(u-T(u)))-f(u,0,0) + f(u,0,0)\ 2 du) 



+2E / \g(u,x(u),x(u - T(u)))\ 2 du. 



By (USD, W> - W - CT and LemmaESl we have 
E(|e( S ,i)| 2 / {pR > 6} ) 

< W\{t-s) J (e\x(u)\ 2 +E\x(u-T(u))\ 2 ^Jdu + 4(t- s)E \f(u,0,0)\ 2 du 
+2Ci y (i + E|x(u)| 2 + E|x(it - t(u))| 2 ) du < C 4 (< - s), 
where C4 is independent of s and i. A combination of this expression and (|3.14p leads to 

E(\e(s, t)\ 2 ) = C A (t - s) + ^ + ^=^. (3-15) 

Therefore, for any given e > 0, we can choose 5 and R such that 

2P +1 6A 1 (p-2)A 1 . 
^^ < 3 6 ' ^/(p-a)^ < 3 6 ' (3 - 16) 

and then choose t — s sufficiently small such that C±(t — s) < |e. Hence, we have 

lim sup E(\x(t) - x(s)\ 2 ) = lim sup E(|e(s, t)| 2 ) = 0. 

t ^ s ^°a<s<t<b t ^ s ^°a<s<t<b 

The proof is complete. 

Lemma 3.8 Assume that problem \2. 1\) € ST>(a, f3, 71, 72). Then for any t\, t 2 ■ a < t\ < t 2 < 

b, 



nx(t 2 ) - y(t 2 )\ 2 < e^C^W+^+^^W+^C*))*!^^) _ y ( tl )\ 2 

(/3(s) + 7i(s)7 2 (s) + 72( s )) e / s t2 ( 2a («)+/ 3 («)+^(«hiH+Ti 2 («))^ds (3.17) 



*2 
fl 



sup E|x(0) -y(0)\ 2 , 

4 0) (ti,t2)<e<i 2 -^i 0) 

where y(t) is a solution of any given perturbed problem 113. fy) . 
Proof. Let 

V{t,x{t)) =p(t)(\x(t)\ 2 +5q(t)), h<t<t 2 , (3.18) 

where 

_ e /*(-2a(«)-,8(«)-7i(«)72(u)-7f(«))dM ) 
q(t) = -(pit))' 1 f (p(u) + 7l(«)7a(«) + 72 2 (^))p(^ 



(3.19) 



(5 is a constant to be determined. For convenience, from now on we use o~(t) and g(t) to denote 
2a(i) + /3(t) + 7i(072(0 + 7i(*) and + 7l(*)7a(t) + 71(0, respectively, that is 



a(t) = 2a(t) +0(t) + Ji(t) 72 (t) + j 2 (t), 
(t) = (3(t)+ 7l (t) l2 (t) + 7 2 (t). 



2 U \ ( 3 - 2 °) 



Then we have p'(t) = -a(t)p(t), (p(t)q(t))' = -p(t)g(t). By (pn8j) . ([37L9l) and the Ito formula, 
one can derive that, for a < t\ < t 2 < b, 

EV(t 2 , x(t 2 ) - y(t 2 )) = EV(h, x(h) - y(h)) 
+ It 2 { ~ a ® p ® E \ x ® ~ y(t)|2 " 

+ 2p(t)E^(x(t) - y(t), f(t, x(t),x(t - r(t))) - f(t, y(t), y(t - r(t)))> 

+p{t)E\g{t, x{t),x(t - r{t))) - g{t, y(t),y(t - r{t)))\ 2 ^dt 

< EV(ti,ar(ti) - y(h)) + J * { - a(t)p(t)E\x(t) - y(t)\ 2 - Sp(t)g(t) 
+2p(t)EK(x(t) - y(t), f(t, x(t),x(t - r(t))) - f(t, y(t),x(t - r(t)))> 
+ 2p(t)ESR(x(t) - y(t), /(t, y(t),x(t - r(t))) - /(t, y(t), y(t - r(t)))> 

+p(t)E| 5 (t, x(t), z(t - r(t))) - g(t, y(t), y(t - T{t)))\ 2 )dt. 



Using the conditions (|2.3p - (|2.5p . we obtain 

EV(t 2 , x(t 2 ) - y(t 2 )) < EV(t!,x(ti) - y(h)) 
+ [ 2 {- <?(t)p(t)E\x(t) ~ y{t)? ~ Sp(t)g(t) 

+2p(t)a(t)E\x(t) - y(t)\ 2 + 2p(t)P(t)E(\x(t) - y(t)\\x(t - r(t)) - y(t - r(t))|) 



+ p(t)E( 11 (t)\x(t) - y(t)\ + l2 (t)\x(t - r(t)) - y(t - r(t))| 



dt 



t-2 



< EV(t 1} x(h) - y(h)) + jf | - <7(t)p(t)E|x(t) - y(t)\ 2 - Sp(t)g(t) 

|2 I „/ -l\ Q ( ,l\ I TV? I ™ / J.\ _ , / \ I 2 I TTJi I / < _/ + \ \ „ + / J. \ \ I 2 



+2p(t)a(t)E|x(t) - y(*)r + p(t)/3(i)(E|x(i) - y(t)\ 2 + E\x{t - r(t)) - y(t - r(t))|' 
+ p(t) 1 j(t)E\x(t) - y(t)\ 2 +p(t) 1 2 {t)E\x(t - r(t)) - y(t - r(t))\ 2 

+2p(t) 11 (t) l2 (t)E(\x(t) - y(t)\\x(t - r(t)) - y(t - r(t))|) left 



which yields 



EV{t 2 ,x(t 2 )-y(t 2 )) 

< EV(t u x(h) - yih)) + p(t)l3(t)(E\x(t - r(t)) - y(t - r(t))\ 2 - 8)dt 

+ J? P{t)l 2 2 (t) (E\x(t - r{t)) - y(t - r(t))\ 2 - s) dt 

+ J? P(*)7l (*)-*(*) (E|x(t - r(t)) - y(t - r(t))| 2 - <5)cft. 



(3.21) 



10 



Lemma l331 implies that E sup |x(i)| 2 < +oo,E sup |y(t)| 2 < +00. Consequently, sup E|: 

a—r<t<b a~r<t<b a~r<t<b 

y(t)\ 2 < +00. Let 5 = sup E\x(t) - y(t)\ 2 . It follows from (ET2T1) that 

/4<V,t 2 )<t<t 2 - A 4 0) 

EV(t 2 ,x(t 2 ) - y(t 2 )) < EV(h, x(h) - y(h)) . (3.22) 
The required estimate (|3.17p now follows from (|3.22[) . 
Lemma 3.9 Assume that problem h2.1\) € ST>{a, ft, 71, 72). Then 

lim sup E\x(t) -x n (t)\ 2 = 0, (3.23) 

n ^°°a-T<t<b 

where x n (t) is a solution of the ^-perturbed problem {2.8\) . 

Proof. For any given natural number n, we can choose a natural number q sufficiently large 
such that /x := (b — a) j q < -. Let 

t\ = a + (i — 1)//, t 2 = a + ip,, i = 1, 2, • • • , q, 

ao = max < max <r(i), Of, /?o = max g(t), 70 = max < max (<r(i) + £>(£)), 1 f 

I a<t<b J a<t<b L a<t<b ) 

V{t,x{t)) = p(t)(\x(t)\ 2 + 5q(t)), h<t<t 2 , 

where 

<5 = e n + sup E|x(6») -x n (6>)| 2 , e„ = sup ( sup E\x(9) - x(t - -)| 2 ), (3.24) 

a-T<e<h a<t<b t-±<e<t H 

n — — 

p{t) and q(t) are defined by (|3.19p . a(t) and g(i) are defined by (|3.20p . For t\ < t < t 2 , we can 
obtain the following estimate in a similar way to the proof of Lemma 13.81 



(3.25) 



EV[t, x(t) - x n (t)) < EV{t x ,x(t x ) - x n {t x )) 
+ H p{8)P(s)(e\x(s - r(s)) - xt t] {s - r(s))\ 2 - S^ds 
+ II P{s)l 2 2 {s) (E\x{s - t{s)) - xt l \s - r(s))\ 2 - 5)ds 
+ ftl p(shi(sh 2 (s)(E\x( S - t(s)) - xt l \s - t{s))\ 2 - 5)ds 
<EV(ti,ar(ii)-x n (ii)) 
+ ( sup Elx(u) — xi n '^(n)| 2 — 5) p(s)g(s)ds. 

Furthermore 

sup E|x(n) — x^ n '^(u)| 2 < sup E\x(u) — x n (u)\ 2 + sup E\x(u) — x n (t )| 2 

a-r<u<t a _ r < u < t _i t _± <u<t n 

— — n n — 

< sup E\x(u) — x n (u)\ 2 + sup ( sup E\x(u) — x n (t )\ 2 ^ < 5, 

a-T<u<t-- a<t<b ^ t ~^-<u<t U 

— — n n — — 

where 5 is defined by (|3.24|) . It follows from (|3.25|) that 

EV(t, x(t) - x n (t)) < KV(ti,x(h) - x n (ti)) , 
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that is, 

E\x(t) - x n {t)\ 2 < e 1 *! ff{u)dtt E\x(ti) - x n (h)\ 2 

+ f g{u)e^ s)ds du(e n + sup E\x{9) - x n {9)\ 2 ) 

Jti ^ a-T<9<t! ' 

< (jH^ u)du + f e {u)e^^ s)ds dn] sup E\x{9) - x n {6)\ 2 + ( I e {u)e& a ^ ds du) 



a-T<e<t\ v Jti 

1 + f (a(u) + Q{u))ef> (s)ds du) sup E\x(9) - x n {9)\ 2 + ( f e {u)e^ s)ds du) i 
< (l + 7o ^ b -^) sup E\x(9)-x n (9)\ 2 + f3 0f ie a °^e n 

a—T<6<a+(i~l)fi 

for all t € [a + {i — 1)^, a + i//], i = 1, 2, • • • , q. Consequently, 
sup E\x(9) -x n (9)\ 2 

a—T<8<a+i(j, 

= max{ sup E\x(9) - x n (6)\ 2 , sup E\x(9) - x n (9)\ 2 } 

a-T<8<a+(i-l)n a+(i-l)^<9<a+ifi * 

< (l + 7 oMe ao{b - a) ) sup E\x(9) - x n (9)\ 2 + f3 0f ie a °( b - a h n 

a—T<9<a+(i—l)fi 

for % = 1, 2, • • • , Therefore, 

sup E|z((9) -x n (6>)i 2 = sup E|z(0) - x n (9)\ 2 

a—T<8<b a—T<8<a+qp, 



<C« sup E\x(9) - x n {9)\ 2 + C * , 1 J 3 ^ b - a h n (3.26) 



= ^(C*-l)e n , 
7o 

where = 1 + 7o//e a °( b_a ). By Lemma 1X71 we have e n = sup sup E|z(0) — x(t — ^)| 2 —?■ 

a<t<bt-l-<8<t 

n — — 

0, as n — >■ oo. Let n — > oo and take into account that 

= ]im L + 70(l>-a) e "'^y = eT „,^,-., 

n— >oo ^ n— >oo \ Q I 

Then (|3.26p leads to the relation (|3.23p . The proof is complete. 

Theorem 3.10 Assume that problem [2A\) G SV(a, (3, 71, 72). lei c = max (2a(t) + 2(3(t) + 

a<t<6 

7l 2 (t) + 2 7l (t) 72 (t) + 7 2 (t)) . Then Vt € [a, 6], 

E\x(t) -y(t)\ 2 < e c (*~ a ) sup E|f(0) - ??(#)| 2 , i/c>0, (3.27) 

a—r<9<a 

E\x(t) - y(t)\ 2 < sup E|f(0) -??(#)| 2 , ifc<0, (3.28) 

a— r<9<a 

where y(t) is a solution of any given perturbed problem 
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(|3.27p and (|3.28p mean that problem (|2.ip is stable in mean square and contractive in mean 
square, respectively. 

Proof. We divide the proof into two cases: [Jf* > and /x^ = 0. 

Case A: fxf 1 ' > 0. In this case, we can obtain the desired result in a similar manner 
as in the proof of Theorem 2.1 in [8]. In fact, replacing a(t), f3(t), \\y(t) — z(t)\\ in [8] with 
2a(t) + 0(t) + 7l(<)72(t) + 7i(t),P(t) + 7l(<)72(*) + 7 2 2 (0^|x(t) - y(t)\ 2 , respectively, using 
Lemma 13.81 and following the proof of Theorem 2.1 in [8], we can obtain either (|3.27[) or ()3.28p 
immediately. 

Case B: /if ) = 0. Note that 

E\x(t) - y(t)\ 2 < E{\x{t)-x n {t)\ + \x n {t)-y n {t)\ + \y(t)-y n (t)\) 2 

= E\x(t) - x n (t)\ 2 + E\x n (t) - y n (t)\ 2 + E\y(t) - y n (t)\ 2 

+2(E(\x(t) - x n (t)\\x n (t) - y n (t)\) + E(\x(t) - x n (t)\\y(t) - y n (t)\) 

+E(\x n (t)-y n (t)\\y(t)-y n (t)\)) 
< E\x(t) - x n (t)\ 2 + E|x n (t) - y n {t)\ 2 + E\y(t) - y n (t)\ 2 



+2((E\x(t) - x n (t)\ 2 y /z (E\x n (t) - y n (t)\ 2 ) 

+ {E\x(t) - x n (t)\ 2 ) 1/2 (E\y(t) - y n (t)\ 2 ) 1/2 
+ {E\x n (t) - y n (t)\ 2 ) 1/2 (E\y(t) - y n (t)\ 2 ) 1/2 ) , 



where y n (t) is the solution of the —perturbed problem 



dy(t) = f n (t,y(t),y(-))dt + g n (t,y(t),y(-))dw(t), t€ [a,b], 
y(t) = r)(t), te[a-T,a], 



(3.29) 



of the problem ()3.2p . f n ,g n are defined by (|2.6j) . It is known that, problem (|2.ip S ST>(a, f3, 71 , 72) 
implies that problem (j2~lj) G 5X>(a, /?, 71, 72). It follows from §ZM) that /2^ 0) > ± > 0. There- 
fore, by case A, for E\x n (t) - y n (t)\ 2 , either (13T2T)) holds if c> or QSjggp holds if c < 0. Letting 
n — > +00 and using Lemma l3.9| we can obtain the desired estimate of E|x(t) — y(t)\ 2 in this 
case. 



Corollary 3.11 Under the assumptions of Theorem \3.1(A Suppose f(t,0, 0) = and g(t, 0, 0) = 
0, then Vt € [a, 6], 

E\x(t)\ 2 < e c (*- a ) sup E|£(0)j 2 , i/c>0, 

a— r<0<a 

E\x(t)\ 2 < sup E|£(#)| 2 , i/c<0. 

a— r<0<a 

Lemma 3.12 Suppose problem \2. 1\) € 5D(a, /3, 71, 72), and i/iai 
2a{t) + p{t) + Ii(t)l2{t) + ll{t) <a <0, 

/3(f)+7iW72(f)+7|W < „ / 1 V * € [°» 6 J' ( 3 - 30 ) 

|2a(i)+0(t)+7i(t)72(*)+7?(*)l ~ ' 
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where ao and v are constants. Then for any given ci, C2, C3 : a < c\ < C2 < C3 < 6, we /lave 

E|z(t) — 2/(i)| 2 < (1/ + (1 - ^)e Q °( C2 - C1 )) sup E|x(0) - y(fl)| 2 , 

V 7 4 0) ( C i,c 3 )<e<c 2 (3.31) 

Vt G [c 2 ,c 3 ]. 

Proof. We divide the proof into two cases: /Xj > and xxf^ = 0. 

Case A: > 0. In this case, replacing a(t),/3(t), \\y(t) - z(t)\\ in [8] with 2a(t) + /3(i) + 
7l(*)72(*) + 7i(*)> + 7l(*)72(*) + T|(*),E|x(t) - y(t)| 2 , respectively, using Lemma [321 and 
following the proof of Lemma 2.3 in [8], we can obtain (|3,3ip immediately. 

Case B: /x^ = 0. In this case, we can obtain the estimate (|3.3ip in a similar manner as in 
the proof of Case B of Theorem 13.101 

3.2 Infinite interval 

Let us now proceed to discuss the equation (|2.1|) which satisfies conditions (|2.2|) - (|2,5|) but the 
integration interval [a, b] replaced by [a, +00). Accordingly, interval [a — r,b] is replaced by 
[a — t, +00). We introduce the symbol ST>(a, (3, 71, 72) to denote this class of problems. 

Theorem 3.13 Assume that problem h2.1\) G SD(a, /?, 71,72), end 

lim (t - r(t)) = +00, sup (2a(t) + f3(t) + 71 (072(0 + 7i 2 (0) = "o < 0, 

<^+°° a<t<+oo 

/3(t)+7i(072(t) + 7 2 2 (0 =y<1 
a <t< +00 \2a(t) + P{t) + 71 (t) l2 (t) + 7? (01 

Then, for any given constant /x > 0, there exists a strictly increased sequence {tk} which diverges 
to +00 as k — > +00, where to = a, such that 

sup E\x(t) - y(t)\ 2 < C* +1 sup E\£(t)- V (t)\ 2 , k = 0, 1, 2, • • • , (3.32) 

ife^t<itfc-l-i a— r<t<a 

where C lt = v + (1- v)e a °^ G (0, 1). Further, 

lim E\x(t) -y(t) | 2 = 0. (3.33) 

t— >+oo 

Proof. It is obvious that (|3.32p implies (|3,33p . So, only the proof of ()3.32p is required. First 
we construct a sequence {t^} by induction. Let to = a. Suppose that tj. is chosen appropriately, 
where k > 0. Because lim (t — r(t)) = +00, there exists a M such that for all t > M, we have 

t— >+oo 

i — T(i) > ifc and therefore /x 2 °^ (M, +00) > tfc. So we can choose ifc + i = M + /x and have the 
relation 

tjfc < (4 ° ~ A*? +°°) < - A* < (3-34) 
Using (|3.34p and Lemma f3. 12 1 we get 

sup E\x(t) - y(t)\ 2 < {v + {l-v)e aolx ) sup E\x(t) - y(t)\ 2 

tk<t<t k+1 $\t k -^t k+1 )<t<t k 

< C M sup E\x(t) - y(t)\ 2 < . . . < C* +1 sup E|£(t) - r?(t)| 2 . 

tfc_l<t<tj; a—r<t<a 

The proof is complete. 
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Corollary 3.14 Under the same conditions as Theorem \3.13[ Furthermore, suppose that f(t, 0, 0) = 
0,s(t,0,0) = 0, then 

sup E\x(t)\ 2 < C* +1 sup E\£(t)\ 2 , fe = 0, 1,2, • • • , 

tk<t<tk+l a—T<t<a 

lim E|x(t)| 2 = 0. 

t— >+oo 

Remark 3.15 J2[/ discussed the stability of nonlinear stiff Volterra functional differential 
equations in Banach spaces. Theorem \3.10\ and Theorem \3.13\ can be regarded generalizations 
of Theorem 2.1 and Theorem 2.2 of /#/ restricted in finite- dimensional Hilbert spaces C d and 
finitely many delays to the stochastic version, respectively. It should be pointed out that function 
f is required to be locally Lipschitz continuous in this paper whereas the corresponding condition 
is not required in J2[/. It is known that local Lipschitz continuity is not a strong restriction. 



(3.35) 
(3.36) 



3.3 Examples 

System (|2.ip includes the following three classes of SDDEs as special cases 

• SDDEs with constant delays: r(t) = r. 

• Stochastic pantograph equations: t — r(t) = qt, where < q < 1 is a constant. 

• SDDEs with piecewise constant arguments: t—r(t) = [t — i\, where [t\ denotes the largest 
integer number less than or equal to t, i is a nonnegative integer. 

Therefore, Theorem 13.101 Theorem 13.131 Corollary 13.111 and Corollary 13.141 are valid for the 
three classes of SDDEs mentioned above. 



Example 3.16 Consider the following linear SDDEs 

dx(t) = (^i(t)x(t) + A 2 (t)x{t - r(t)) + F(t))dt 

+ (B 1 (t)x(t) + B 2 {t)x{t - r(t)) + G(t))dw(t), { > 

where Ai(t) , A 2 (t) , Bi(t) , B 2 (t) G C dxd ,F(t),G(t) G C d are continuous with respect to t, w(t) is 
an 1-dimensional Wiener process. For the problems {3.31 ), it is easy to verify the conditions 

(KM) -(KM) with 

a(t) = A^r^ lW , f3(t) = \A 2 (t)\, 7 i(i) = 7 2 (i) = \B 2 (t)\, 

A^(t)+A 1 (t) 

where A max 2 denotes the largest eigenvalue of the Hermite matrix diW+^iW _ 

Applying Theorem 13.101 and Corollary 13.111 to (j3.37j) , we have 
Corollary 3.17 The solutions of {3.37]) satisfy 

E\x(t) - y(t)\ 2 < e c ^ sup E|£(0) - r](6)\ 2 , t G [a,b], if c > 0, 

a—r<9<a 

E\x(t) - y(t)\ 2 < sup E|£(0)-J7(0)| 2 , t G [a, b], if c < 0, 

a~T<8<a 
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where x(t),y(t) are the solutions of \3.31ty corresponding to the initial functions £(t) and rj(t), 
respectively, 

+ 2|A 2 (t)| + (|B 1 (t)| + |B 2 (t)|)' 



nriax 



a<t<b 

Furthermore, if F(t) = and G(t) = 0, then the solutions of {3.31^ satisfy 
E\x(t)\ 2 < e^-^ sup E|£(#)| 2 , t G [a, b] , if c> 0, 

a—T<6<a 

E\x(t)\ 2 < sup 3E|^(6») [ 2 , t e [a,b], if c < 0. 

a—T<8<a 

Applying Theorem 13.131 and Corollary 13.141 to (|3.37|) , we have 
Corollary 3.18 If lim (t - r(t)) = +oo, 

t— >+oo 

sup 2A max 2 + \A 2 (t)\ + |fli(i)||B 2 (t)| + l^^)! 2 < 0, 



a<t<+oo 



„ m |A 2 (f)| + |i? 1 (f)||B 2 (t)| + |B 2 (f)p 1 

blip A*(t)+A!(t) v 



a<t<+oo | , 2" 



2A max 2 +1^2(4)1 + 151 (t)||B 2 (t)| + |Bi(t) 



- 



then the solutions of ^3.37\ ) satisfy lim E[a:(t) — y(i)| 2 = 0- Furthermore, if F(t) = and 

t— i> + 00 

G(t) = 0, i/ien i/ie solutions of fOTU safc/y lim E|x(t)| 2 = 0. 

t— >+oo 

In particularly, if d = 1, ^4i(t), ^4. 2 (t), Bi(t), B 2 (t) are constants, that is, vli(i) = Ai,A 2 (t) = 
A 2 ,Bi(t) = B\,B 2 {t) = B 2 , then the problems (|3.37p are mean-square stable if 

^A 1 + \A 2 \ + ^\B 1 \ + \B 2 \y <0. (3.38) 

Remark 3.19 For linear scalar SDDE with constant delay and linear scalar stochastic panto- 
graph equation, fTl\ [7Tj) / and J2C obtain the condition 13. 38]) . respectively. Therefore, specializing 
Corollary \3.14\ to the case of linear scalar SDDE with constant delay and linear scalar stochastic 
pantograph equation, respectively, the corresponding result is accordant with that in the literature. 

Example 3.20 Consider the following nonlinear equation 

dx(t) = (Ai(t)x + A 2 (t)x 3 + A 3 (t)Jx 2 (t - r(t)) + 1 + F(tj) dt 

/ \ (3.39) 

+ [Bi(t)smx(t) + B 2 (t)arctanx(t - T(t)) + G(t))dw{t), 



where Ai(t),A 2 (t),As(t),Bi(t),B 2 (t),F(t),G(t) are continuous real-valued functions in t and 
A 2 {t) < 0. It is easy to verify that $3.39\) satisfies the conditions i2.2\) - [275\) with 

a(t) = A 1 (t), P(t) = \A 3 (t)\, 7l (t) = IfliCt)!, l2 {t) = \B 2 {t)\. 

Applying Theorem l3.10l Corollarv l3.11l Theorem l3.13l and Corollarv l3.14l to ()3.39p . we can obtain 
properties of solutions of (|3.39p which are similar to those of Corollary 13.171 and Corollary 13.181 
For the sake of brevity, we do not present the corresponding results for ()3.39p . 

Remark 3.21 The drift coefficient f(t, x, y) of \3. 39\) satisfies local Lipschitz condition and one- 
sided Lipschitz condition with respect to x but global Lipschitz condition. The stability analysis in 
this paper is based on the local Lipschitz condition and the one-sided Lipschitz condition instead 
of a more restrictive global Lipschitz condition. 
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4 Stability of backward Euler method 

In this section, we investigate whether numerical methods can reproduce the contractivity in 
mean square. For the deterministic differential equations, it is known that the contractivity of 
numerical methods is too strong [2J [TJ. The existing theory Q2] shows that only backward 
Euler method and two-stage Lobatto IIIC Runge-Kutta method can preserve the contractivity 
of nonlinear delay differential equations. Therefore, in the stochastic setting, we only focus on 
the backward Euler method instead of other methods. 
For simplicity, from now on, we assume that 

a(t) = a, p(t) = f3, 7i (t) = 71, 72(f) = 72, t€ [0,6]. 

On a finite time interval [a, b], a uniformly partition is defined by 

b — a 

U = a + in, 1 = 0, 1, . . . , h = — . 
The backward Euler (BE) method applied to (|2.ip yields 

^ X n+ \ = X n + hf(t n+ i, X n+ \, X h (t n+ i — r(t n +i))) , . 

+g(tn, X n , X h {t n - T (t n )))Aw n , n = 0, 1, . . . , N - 1, 1 ' a) 

^X h (t) = 7r h (t,C,X 1 ,X 2 ,...,X n ), a-r<t<t n , (4.1b) 

where 7r h is an appropriate interpolation operator which approximates to the exact solution 
x(t) on the interval [a — T,b], X n is an approximation to the exact solution x(t n ), Aw n = 
w{t n+ \) — w(t n ). It is well known that the backward Euler method is convergent with strong 
order only 1/2 for stochastic differential equations. We can use the following piecewise linear 
interpolation 

vh( + \ _ J - t)Xi + (t- ti)X i+1 },U < t < t i+1 ,i = 0, 1, 2, . . . , N - 1, 

[) ~\ Z(t), a-T<t<a. [ ^ Z) 



Applying BE to the perturbed problem (|3.2p we can obtain the corresponding scheme 

Y n+1 = Y n + hf(t n+1 ,Y n+1 ,Y h (t n+1 -T(t n+1 ))) 

+g(t n ,Y n ,Y h (t n -T(t n )))Aw n , n = 0,l,...,N-l, (4.3) 
Y h (t) = ir h (t, T), Y X ,Y 2 , . . . , Y n ), a-T<t<t n . 



For simplicity, for any given nonnegative integer n, we write 

P n = X n -Y n , Q n =max{max E|Pi| 2 , sup E\£(t) - V (t)\ 2 }, n > 1, 

l<i<n a-r<t<a 

Qo= sup El^t)-^)! 2 . 

a—T<t<a 



(4.4) 



Moreover, for convenience, we introduce notations to denote the values of drift and diffusion 
coefficients at specific points. 

f xx (n + 1) = f(t n+1 , X n+1 , X h (t n+1 - r(t n+1 ))), 

fW(n + 1) = f(t n+1 ,Y n+1 ,Y h (t n+1 - r(Wi))), ,s 

fv»(n + l) = f(t n+l ,Y n+l ,X h (t n+1 - T (t n+1 ))), { - 0) 
g xx (n) = g(t n ,X n , X h (t n - r(t n ))), g v»(n) = g(t n , Y n , Y h (t n - r(t n ))). 
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Lemma 4.1 Under the conditions \2. 3\) and &2.4\ )> if {a + f3)h < 1, the implicit equation f^, la ) 
admits a unique solution. 



Proof. Let f(z) := f(-,z,z h (-)), then implicit equation (|4.1ap can be rewritten as 

z = hf(z) + b = hf(-,z,z h (-)) + b, (4.6) 

where z is unknown whereas b and h are known. Inserting the interpolation operator (|4.2p into 
(|4.6|) . we have 

z = hf(z) + b = hf(; z, Iz + b ) + 6, (4.7) 
where < / < 1, Z and 6o are also known. It follows from (j2.3f) . (j2.4|) and ()4.7p that 

R(zi - Z 2 , f(zi) - f(z 2 )) = U{Zt - Z 2 , /(", Zx, lZl + 6o) - /(') 32, ^2 + b )) 

= SR^i - z 2 ,f(-,zi,lzx + b ) - f(-,z 2 ,lzx + b )} 

+0?(^i - z 2 , f(-,z 2 ,lzx + 6 ) - f(-,z 2 ,lz 2 + 6 )) < «ki - ^2 1 2 + /3|zi - 2 2 | 2 . 

The assertion follows immediately from [H Theorem 5.6.1]. 

Theorem 4.2 Assume that problem 112. 1\) G ST>(ct, ,0,71,72)- iei -pf n } On} ^ e se " 
quences of numerical solutions obtained by the backward Euler schemes j^. i| ) and |J3]), respec- 
tively. Write c = 2a + 2/3 + 7^ + 27172 + 7^. 

(%) If c> 0, for any given cq € (0, 1), f/ien we have for he < Co 

E|X n -y n | 2 < e^"- a ) sup E|£(t)-??(t)| 2 , n = 1, 2, • • • , N, (4.8) 

a—T<t<a 

wftere c- T ,ci - max| x-2ha-2hp ' i-2ha-hp > > L - 

(ii) If c < 0, i/ien we have for any h > 

E|X n -y n | 2 < max E|£(t) - rj(t)\ 2 , n = 1, 2, • • • , N. (4.9) 

a—T<t<a 

Note that the stability in mean square (|4.8p and the contractivity in mean square (j4.9|) can 
be regarded as numerical analogs of (|3.27p and (|3.28p for the analytical solution of the problem 
(|2.ip . respectively. 

Proof, (i) By g4]) and (|3~5|> . we have 

Pn+i - Kf> + 1) - / TO (n + 1)) = P n + (<f » - 9 »(n))A %) 

which yields 

|P n+1 | 2 - 2hM{P n+1 ,f xx (n + 1) - /«*(„ + 1)> + ^ir^n + 1) - /««(„ + l)| 2 
= |P n | 2 + 2^<P n , (^(n) - g yy(n))Aw n ) + \ (<f » - <^(n)) A^ n | 2 . 
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Taking expectation and using (|2.3p - (|2.5p and (|4.2p . we get 

E|P n+1 | 2 < E|P n | 2 + 2hm(Pn+i, f xx (n + 1) - /W(n + 1)> + hE\g xx (n) - 5 W N| 2 

< E|P n | 2 + 2hm(P n+1 J xx (n + 1) - /»*(n + 1)> 
+2/iE»(P n+ i, /^(n + 1) - py(n + 1)> + ZiE^^n) - g yy {n)\ 2 

< E|P n | 2 + 2/iaE|P n+ i| 2 

+2hf3E(\P n+1 \\X h (t n+1 - r(t n+1 )) - Y h (t n+1 - r(t n+1 ))\) 
+hE( 7l \P n \ + l2 \X h (t n - T(t n )) - Y h (t n - r(t n ))\) 2 

< E|P n | 2 + 2/ l aE|P n+ i| 2 

+/i/3E|P n+1 j 2 + h/3E\X h (t n+ i - r{t n+1 )) - Y h (t n+l - r(t n+1 ))\ 2 

+Kll + 7 i72)E|P n | 2 + h^l + lll2 )E\X h {t n - r(t n )) - Y h (t n - r(i n ))| 2 

< E|P n | 2 + 2/ l aE|P„+i| 2 + / l /3E|P n+1 | 2 

+h/3 max { max E|P| 2 , sup E|£(i) - v{t)\ 2 } 

l<i<n+l a-r<t<a 

+Mt? + 7i72)E|P n | 2 + h(d\ + 7i 72 )Qn, 

(4.10) 

where we used the piecewise linear interpolation (|4.2p and the following inequality 

E|(l-5)P + 5P+i| 2 <max{E|P| 2 ,E|P +1 | 2 }, < S < 1. (4.11) 
It is clear from (|4.1Up that 

(1 - 2ha - hp)E\P n+1 \ 2 < (1 + hj 2 + /i 7l72 )E|P n | 2 



+/i/3max{ max E|P| 2 , sup E|£(t) - r]{t)\ 2 } + ft( 7 | + 7 i 72 )Q r 

l<i<n+l a-r<t<a 



(4.12) 



We now consider two cases: 



(a) max{ max E|P| 2 , sup E\£(t) - r](t)\ 2 } = E\P n+1 \ 2 , 

l<i<n+l a-r<t<a 

(b) max{ max E|P| 2 , sup E[f(t) - r](t)\ 2 } / E|P n+ i| 2 . 

l<i<n+l a-r<t<a 



In the case of (a), it follows from (|4.12|) that 

(1 - 2ha - 2/i/3)E|P n+1 | 2 < (1 + hrfi + 2h 7x72 + hj 2 )Q n , (4.13) 

which yields 

E|P n+1 | 2 < 1 + + L ^ 7172 + H Q n < c lQn . (4.14) 

In the case of (b), (|4,12p implies that 

(1 - 2ha - h{3)E\P n+1 \ 2 < (1 + /i 7l 2 + /i 7l72 )E|P n | 2 + hf3Q n + h( 7 2 + 7 i 72 )Q n 
< (l + / i /3 + /i 7l 2 + 2/ i7 i 72 + /i 72 2 )Q n , 

which yields 

m w i 1 + v+ , 1„ + 2h lT + hll Qn < ■»<*.. (4.i5) 

1 — lha — hp 
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To summarize, both in the cases we have shown that E|P n+ i| 2 < c\Q n , which yields 

Qn < Qn-1 + E|P n | 2 < (1 + Ci)Q n _i. (4.16) 

By induction, we further obtain 

E\X n -Y n \ 2 = E|P n [ 2 <Q n < (l + ci)Q„_i < ••• < (l + Cl ) n Q 
< e Cin Q = e^"- a) sup E|£(t) - r](t)\ 2 . 

a—T<t<a 

(ii) When c < 0, noting that (jlTHD . (|4"7L5D and 

1 + fry 2 + 2^7172 -f/rg < 1 + + h-fl + 2/t7i7 2 + hj% < 
l-2ha-2h(3 ~ ' l-2ha-h/3 ~ ' 

we have for any h > 

E|X n -y n | 2 <Q n _ 1 <Q n _ 2 <---<Qo= sup E|£(t)-n(£)| 2 . (4.17) 

a— r<i<a 

Therefore we have completed the proof of the theorem. 

Corollary 4.3 Assume that problem 12. 1\) G 5D(a, ,$,71,72). Lei 6e a sequence of nu- 

merical solutions obtained by the backward Euler method fl^. i[ ). Furthermore, if f(t, 0,0) = 
0,s(t,0,0) = 0, tfien 

(%) /or i/ie case 0/ c > ; /or any giwen Co € (0, 1), we aai>e /or /ic < Co 

E\X n \ 2 < e £{tn ~ a) max E|£(i)| 2 , n = 1, 2, • • • , N; (4.18) 

a— r<t<a 

(ii) for the case of c < 0, we have for any h > 

E|X„| 2 < max E|£(t)| 2 , n = 1, 2, • • • , N. (4.19) 

a—r<t<a 

Theorem 4.4 Assume that problem \2. 1\) G 5P(a, /?, 71, 72), and 

lim (t - r(t)) = +00, c = 2a + 2/3 + 7l 2 + 2 7l72 + 7! < 0. (4.20) 

t— >+oo 

Lei {X ra } and {1^} be two sequences of numerical solutions obtained by the backward Euler 
schemes fl^. i[ ) and fl^.3[ ). Taen, 

taere exists a strictly increased positive integer sequence {n^} which diverges to +00 as 
k — > +00, where no = 0, such that for any given h > 0, 

max E\Xi-Yi\ 2 <c% +l sup E[f(t) - r/(t)| 2 , A; = 0, 1, 2, ... , (4.21) 

n fe <i<ra fc+1 a-r<t<a 

waere c 2 - max| i_ 2 / lQ -2/ l /3 ' i-2ha-hp / < x > 
(mJ /or any given h > 0, 

lim E|X n - Y n | 2 = 0. (4.22) 

71— > + 00 
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Note that the contractivity in mean square (|4.2ip and the asymptotic contractivity in mean 
square (|4.22|) can be regarded as numerical analogs of (|3.32|) and (|3.33p for the analytical solution 
of the problem (|2.1|) . respectively. 

Proof. It is obvious that (|4.21|) implies (|4.22|) . and we only need to prove (|4.2ip . By (|4.20p . we 

have 2a + j3 + 7172 + < 0, (/3 + 7172 + 7 2 )/|2a + P + 7i72 + 7il < 1 and c 2 < 1. 

First, as done in 0HE], we can construct a strictly increased sequence of integers {n^.} which 
diverges to +00 as k — > +00, such that 

t - r{t) > t nk+ i, Vt > t nk+1 , 

where n$ = 0. In fact, suppose that nk(k > 0) has been chosen appropriately. Then there exists 
a constant M > t nk such that for all t > M we have t — r(t) > t nk + h since limt_ ++00 (t — r(t)) = 
+00. If M is a node, we let t nk+1 = M, otherwise there exists natural number j such that 
tj < M < tj+i, then we let rik+i = j + 1 and t nk+1 = tj+\. Thus we obtain the required 
sequence {n^} which satisfies 

to < t\ < ■ ■ ■ < t ni < t ni+ i < ■ ■ ■ < t n2 ■ ■ ■ < t nk < ■ ■ ■ . 

For rif. < rt + 1 < n£. +1 , by the second inequality of (|4.1U|) and conditions (I2.3|) - (l2.5p . we have 
E|P n+1 | 2 < E|P n | 2 + 2/iE^P n+1 ,/^ + l)-/^(n + l)> 

+2hEM(P n+1 jy x (n + 1) - fyy(n + 1)> + hE\g xx (n) - g yy (n)\ 2 

< E\P n \ 2 + 2haE\P n+1 \ 2 + 2hpE{\P n+l \\X h (t n+1 - r(t n+1 )) - Y h (t n+1 - r(t n+ i))|) 

+hE(j 1 \P n \ + l2 \X h (t n - r(t n )) - Y h (t n - r(t n ))|) 2 

< E|P n | 2 + 2haE\P n+l \ 2 + h(3E\P n+1 \ 2 + h(3E\X h (t n+1 - r(t n+1 )) - Y h (t n+l - r{t n+ 
+Kl\ + 7i72)E|P n | 2 + h^ 2 + lxl2 )E\X h {t n - tQQ) - Y h (t n - r(t n ))\ 2 , 



which yields 

(1 - 2ha - hp)E\P n+ i\ 2 

< (1 + h 7 2 + h lll2 )E\P n \ 2 + hf3E\X h (t n+1 - r(t n+1 )) - Y h (t n+1 - r(t n+1 ))\ 2 
+Kll + 7 i72)E|X h (t n - r(t n )) - Y h (t n - r{t n ))\ 2 

< (l + /i 7l 2 + / i 7i7 2 )E|P n | 2 + /i/3 max E|P| 2 + h(j 2 + 7l72 ) max E|P;| 2 , 

nfe_i<i<n+l n^_i<i<n 

where we used the piecewise linear interpolation operator fj4.2j) and the inequality ()4.1ip . We 
now consider the following two cases. 
If max E|P| 2 = E|P n+ i| 2 , we have 

UTID 12 S 1 ± + 2/t 7l72 + hrj W |D|2^ W |d|2 

E[P„+i| < — —z max E|Pj| < c 2 max E|P| . 

1 — ZtlCt — Zllp nj._!<i<n nj ; _ 1 <i<n 

If max E|P| 2 7^ E|P n+ i| 2 , we have 

n k _ 1 <i<n+l 

IB- IP ,2^ 1 + ^/3 + Wl + 2/17172 + hf£ W |p|2^ TH.rp.2 

E|P n+ i| < — max IE, I Pj I < c 2 max E|Pj| . 

1 — 2\lOL — hp nf i _i<i<n nj._i<i<n 
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In both cases, we have obtained the inequality 



E|P n+ i| 2 < c 2 max E|P;| 2 , n k <n+l<n k+1 . (4.23) 

n k _ 1 <i<n 

(|4,23p with n = n k reduces to E|P rafc+ i| 2 < c 2 max E|Pj| 2 . By induction, we have 

«fe-i<«<«fe 

max E\Xi-Yi\ 2 = max E|Pj| 2 < c 2 max E|P;| 2 

n k <i<n k+1 n k <i<n k+1 n k _ 1 <i<n k 



< ■■■<4 +L max E|£(t)-7/(t)|' 

a— r<t<a 



The proof is complete. 



Corollary 4.5 Under the same assumptions of Theorem \4-4\ Let {X n } be a sequence of nu- 
merical solution obtained by the backward Euler method J^,i| ), Furthermore, if f{t, 0,0) = and 
g(t, 0,0) = 0, then, 

(i) there exists a strictly increased positive integer sequence {n^} which diverges to +oo as 
k — > +oo, where no = 0, such that for any given h > 0, 

max E|Xi| 2 <c^ +1 sup E[£(t)[ 2 , k = 0,1,2,.... 

n k <i<n k+1 a-r<t<a 

(ii) for any given h > 0, lim E|X n | 2 = 0. 

n— »+oo 

5 SDDEs with several delays 

Consider the following SDDEs with several delays 

dx(t) = f(t,x{t),x(t-Tx (*)),■•• ,x(t-T r (t)))dt 

+g(t,x(t),x(t-T 1 (t)),--- ,x(t-T r {t)))dw(t), t>a, (5.1) 
x(t) = £(t), t G [a-r,a], 

where Tj(t) > 0, i = 1, 2, • • • ,r and max inf (t — Tj(t)) > a — r. All results given in this paper 

l<i<r t>a 

can be extended easily to the case of several delays. For the sake of brevity, we do not present 
the corresponding results for (|5.ip . 



6 Conclusions and future work 

In this paper, we investigate stability of analytical and numerical solutions of nonlinear SDDEs. 
Sufficient conditions for stability in mean square, contractivity in mean square and asymptotic 
contractivity in mean square to the solutions of nonlinear SDDEs are given, which provide a 
unified theoretical treatment for stochastic differential equations, stochastic delay differential 
equations with constant and variable delays, bounded and unbounded delays . Furthermore, it 
is proved that the backward Euler method can preserve the properties of the underlying system. 
The main results of analytic solution in this paper can be regarded as generalization of those 
in [S] restricted in finite-dimensional Hilbert spaces and finitely many delays to the stochastic 
version. We have encountered problems when we tried to obtain a unified framework for general 
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stochastic functional differential equations. It is worth noting that whether the results in [8] 
can be extended to general stochastic functional differential equations or not. One area for the 
future work is to give a positive or negative answer for the question. Neutral stochastic delay 
differential equation (NSDDE) is more general than stochastic delay differential equation. It 
is interesting to investigate whether the theory of this paper can be extended to NSDDEs and 
corresponding numerical methods. It will also be our future work. 
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